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Abstract 

To understand the mechanism allowing for the long-term storage of 
excess energy in proteins, we study a Hamiltonian system consisting 
of several coupled pendula in partial contact with a heat bath. It is 
found that energy storage is possible when the motion of each pendulum 
switches between oscillatory (vibrational) and rotational (phase-slip) 
modes. The storage time increases almost exponentially to the square 
root of the injected energy. The relevance of our mechanism to protein 
motors is discussed. 

Proteins are among the most important biopolymers for living systems. 
They transform chemical energy to mechanical energy, and vice versa, and 
contribute to biological functions. However, the question of how proteins work 
dynamically remains unanswered. Recently, a noteworthy experiment concern- 
ing protein motors was performed In this experiment, the working process 
of a single molecule was directly investigated. The results suggest that proteins 
often store energy obtained from a reaction with ATP (adenosine triphosphate) 
and use it later (e.g., for enzymatic reactions with other proteins). The inter- 
val for energy storage was found to sometimes be very long, up to the order 
of seconds, while typical timescales for normal vibrations are several picosec- 
onds. How can proteins store excess energy for such a long time, somehow 
overcoming the relaxation process toward thermal equilibrium? In order for a 
protein to store energy for a sufficiently long time, energy must be absorbed 
into a certain part of the protein, in accordance with its own dynamics. Fur- 
thermore, some characteristic type of dynamics is required to store the excess 
energy without losing it to the surrounding aqueous solution. 
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As a first approach to understand the working mechanism of proteins, we 
construct a Hamiltonian system (in partial contact with a heat bath), which 
stores energy for a given time span, in spite of its eventual relaxation to thermal 
equilibrium on a much longer timescale. In this Letter, we adopt a system 
consisting of several coupled pendula 0, |], each of which possesses two modes 
of motion, oscillation and rotation. We clarify the characteristic dynamics 
necessary for energy storage in connection with the coexistence of these two 
modes, and also show that partial contact with the heat bath is necessary for 
this storage. The relationship between the storage time and the injected energy 
is obtained, and is experimentally verifiable. The relevance of our results to 
protein motors is also discussed. 

A protein consists of a folded chain of amino acids that assumes a globular 
shape [0. The main chain is accompanied by side chains arranged around it, 
with each side chain hanging on the main chain similar to a pendulum. Some 
side chains are gathered in a globular shape, and a certain assembly of them 
can play an important role in the function of the protein. As an abstract model 
for the angular motion of side chains in such a functional assembly, we choose 
a system of coupled pendula. In particular, we study the idealized case of N 
identical pendula equally coupled to each other. Here, the oscillation of the 
pendula corresponds to the vibration of the side chains. With this simple model 
as an example, we demonstrate that long-term storage is generally possible 
in a class of Hamiltonian systems Our study is restricted to a classical 
mechanical description, since quantum mechanical effects are believed to be 
irrelevant to protein dynamics (except for the choice of the potential). 

Our Hamiltonian is given by 

N 2 N 

h=k+v=j:^+y: v{9,,9,), (1) 

where pi is the momentum of the i-th pendulum. The potential V is con- 
structed so that each pair of pendula interacts through their phase difference 
with an attractive force to align the phases 0, |^: 

= 2(2^^^ " cos(27r(^, - 6,))}. (2) 
Hence, the evolution equations for the momentum pi and the phase 6i are given 

by 

1 ^ 

= ^T.MMo^-o,)), (3) 

i=l 

e, = p,. (4) 

In this model, all pairs of pendula interact identically. If the protein were a 
straight chain without folding, it could be modeled by a one-dimensional chain 
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of pendula. In reality, due to its globular shape, a global interaction occurs 
between pendula. Although the assumed "mean-field coupling" with uniform 
strength represents an extreme simplification, the resultant model can capture 
some general characteristics of protein dynamics. 

For the (dynamical) function of proteins, it is essential that they exist in 
aqueous solution, which functions as a heat bath. Since the hydrophobic part 
of a protein molecule is segregated from the aqueous solution, contact with the 
heat bath is restricted. Accordingly, we define the equation of motion so that 
only some pendula in the system contact the heat bath: 

PjTT 

P^ = —QQ;^-lPi + ^m {l<Nh). (5) 

Here, the heat bath is described by the Langevin equation, in which T rep- 
resents its temperature, 7 is a relaxation coefficient, and ^j(t) is a Gaussian 
random form satisfying {^iit)) = and (^j(ti)^j(t2)) = 5ij^it2 — ti), with (■) as 
the temporal average. An important point here is that contact with the heat 
bath is restricted to the few pendula satisfying i < N^. This restriction may 
be interpreted as only allowing for interaction of the hydrophilic part with the 
aqueous solution. Note that this restricted heat bath (i.e., partial contact with 
the heat bath) is sufficient for realizing thermal equilibrium. 

Here, we briefiy review the typical behavior of eqs. (^ and (^ in the 
conservative case 0, |^, |^, without coupling to the heat bath. When the 
total energy is small, all pendula vibrate almost synchronously (i.e., \9i — 9j\ 
and \pi — pj\ are small for all i, j). The typical timescale of vibration is ~ 1 
- 10, while there is collective motion giving rise to a periodic change in the 
degree of synchronization, with a timescale of ~ 10^ - 10^. 

As the total energy increases, rotational (phase-slip) behavior of pendula 
begins to appear. For a sufficiently large total energy, almost all pendula ex- 
hibit rotational motion. Here, each pi changes slowly over time, with relatively 
large values of \pi—pj\. In this case, the pendula rotate almost freely, since the 
correlation of their phase with those of other pendula cannot be maintained, 
and the interaction term cancels out even with a short time average, except in 
the rare situation that two momenta take very close values. The timescale of 
rotation is 0(1). 

In the medium energy regime, the rotational motion of a single pendulum 
appears intermittently from an assembly of vibrating pendula (as displayed 
in the upper part of Fig. Once a pendulum starts to rotate, it typically 
continues to rotate over many cycles, which is longer than the typical timescale 
of each pendulum's vibration (and rotation). We emphasize here that the 
effective interaction for a rotating pendulum is much weaker than those for 
vibrating ones, as mentioned above. 

In the thermodynamic limit with N ^ 00 p, ||, 0], there appears a phase 
transition at the energy E^/N = 0.0190 (T^ = 0.0127), where the vibrational 
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mode is dominant in the 'solid'-like phase and the free rotation is dominant in 
the 'gas'-hke phase. If is not sufficiently large (as is the case studied here), 
the transition point is blurred by the finite size effect, and the fluctuations 
are large around the transition point. Hereafter, we refer to this temperature 
range (0.01 ~ T ~ 0.03) as the medium temperature range. 

Now, we consider the Hamiltonian system in contact with the heat bath 
described by eq. (^, where the total energy is now time-dependent. It is 
found that the fluctuations of energy are sometimes correlated with the dy- 
namics of the pendula. For instance, in the medium temperature regime in 
partial contact with the heat bath, the intermittent appearance of rotating 
pendula, similar to the Hamiltonian dynamics, becomes accompanied with a 
large fluctuation of the total energy (see Fig. |^). Here, deviation toward a 
larger total energy is supported by the concentration of energy in one (or few) 
pendulum [g]. 

While local fluctuations are thus dependent on the dynamics, the system 
reaches thermal equilibrium for a sufficiently long timescale, irrespective of the 
values of 7, Nh (> 1), or the temperature. Note that the equihbrium property 
as a canonical ensemble is identical for any number of N^, even for Nh <C 
A^. When Nh = N (full contact), the timescale to reach thermal equilibrium 
is determined completely by the value of 7, independently of the dynamical 
properties of the system. On the other hand, in the case of partial contact, the 
relaxation process is affected by inherent Hamiltonian dynamics. It progresses 
beyond the timescale F"^, where F = 'jN^/N gives the dissipation rate of 
the system. Even in this case, the dynamics of the pendula corresponding to 
i ^ Nh are governed by the timescale 7"^ (< F~^). 

As the next step, we discuss the energetic behavior of the system far from 
equilibrium. Consider a special enzymatic event such as ATP attachment or 
its reaction. With such an event, proteins are moved far from equilibrium. We 
study how such an event is related to the storage of energy. 

Due to such a 'reaction' event, some portion of the protein is forced far 
from the previous equilibrium state. This situation is modeled by the addition 
of an instantaneous kick to a certain pendulum at the reaction event. The 
kicked pendulum comes to possess a larger amount of momentum and kinetic 
energy than other pendula §]. 

An example of temporal evolution after the kick at t = to is shown in Fig. 
H, where the system with A^ = 10 and Nh = 1 is adopted; that is only the first 
pendulum is in contact with the heat bath. Kinetic energy with the amount of 
£^0 is added to the kicked pendulum. In Fig. |^, the high-energy state continues 
up to 2 X 10^ while the typical relaxation time F~^ of the heat bath is equal to 
10^. There, the kicked pendulum continues to rotate in isolation, maintaining 
the large energy, while it is affected only slightly by the other pendula over a 
long time interval. This allows for long-term energy storage. 

In order to study the lifetime of energy storage, we define it as the interval 
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from the kick until the relaxation of the total energy (see Fig. 0). This is 
determined by the time at which the total kinetic energy K decreases to 
NT/2, that is the value at thermal equilibrium. The results, however, do not 
depend on the specific choice of the relaxation time. 

The distribution of the lifetime is shown in Fig. ^, where evolutions for 
200 stochastic processes ^j(t) are sampled. For comparison, we show the dis- 
tribution for two cases with the same value of F: = 1 with 7 = 10~^ and 
Nfi = N{= 10) with 7 = 10~^. It is noted that the distribution of the lifetime 
is quite different for the two cases. For the case of partial contact, Nh = 1, the 
typical lifetime of the energy storage is very long and reaches 10^, in contrast 
with 10^ for the A^";, = case. Moreover, the long-term energy storage is 
obtained as far as the kicked pendulum is not in contact with the heat bath, 
even if ~ as long as 7^ N. This result suggests that partial contact 
with the heat bath is necessary for long-term energy storage. 

Although the thermal equilibrium properties are similar in the above two 
cases {Nfi < N and = N), a large difference appears when the system is 
placed far from equilibrium. The pendula i < N^, interacting directly with 
the heat bath, cannot rotate freely over the I/7 timescale (in contrast with 
the Hamiltonian dynamics) and the dynamics is replaced by Brownian motion 
due to the heat bath. In the case of partial contact, a long duration of rota- 
tion is possible for pendula corresponding to i > A^^, if the pendula possess 
sufficiently large energy to remain far from equilibrium. The pendula there 
become free from the thermal effect, due to the distinctively weaker interac- 
tion with the other pendula, and follow nearly pure Hamiltonian dynamics. 
Although the prototype of this mechanism is observed as a large fluctuation 
around equilibrium in the medium temperature regime (see Fig. |1]), it works 
well far from equilibrium. 

The lifetime of energy storage increases with the increase of the kicked 
energy Eq, injected to a single pendulum. Figure | shows the relationship 
between the average lifetime and Eq. The cases in partial contact = 1 with 
two different 7 values are compared with the Hamiltonian (i.e., microcanonical) 
case without the heat bath |T0[. In the Hamiltonian case, we find the following 
relation 

{tR - to) oc exp(ay^), (6) 

between the average lifetime and Eq with (■) as the ensemble average, and a as 
a temperature-dependent constant. In a dissipative case with the heat bath, 
a similar increase of the relaxation time is maintained, although there exists 
slight suppression of the lifetime at a high energy. 

In Hamiltonian systems, the exponential law of relaxation time is known as 



the 'Boltzmann- Jeans conjecture' []TT], |T2|, |T^, [IJ]. The energy relaxation be- 
tween slow and fast modes generally requires a considerable amount of time. 
The conjecture is confirmed in numerical experiments for a classical gas of di- 



atomic molecules |]T5|, and some analytic estimates for a generalized version of 
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the conjecture are given within a classical perturbation theory of Hamiltonian 



systems |p.6|| . In our case, there is only slight energy exchange between the fast 
rotational and slow vibrational modes. In the above diatomic molecules with 
translation, rotation, and vibration, when it is sufficiently high, the rotational 
energy almost freezes, and the transfer of the energy to the translational mode 
requires a time exponential to the angular momentum, similar to eq. (P). We 
may expect that this exponential form is generally valid, if the excited mode 
has a weaker interaction with other modes as the excited energy becomes 
higher. 

Thus far, the exponential law has been proven for a class of Hamiltonian 
systems. Although our model has partial contact with the heat bath, a similar 
increase of the relaxation time is confirmed. Moreover, the same increase is 
obtained for a system with the same F value, for any value of (up to iV — 1), 
as long as the rotating pendulum does not come into direct contact with the 
heat bath. On the other hand, if all pendula are in contact with the heat bath, 
the energy of the rotating pendulum decays much faster than the case with 
partial contact depicted in Fig. |. The typical lifetime for this full contact 
case is of the order F~^ over a wide range of £"0. The system in partial contact 
with the heat bath retains a mechanism for slow relaxation, similar to the 
Hamiltonian system. 

What type of Hamiltonian is required for long-term energy storage? Con- 
sider a system with two possible phases (e.g., gas-like and solid-like phases). 
Here, some elements possessing higher energy "melt" and come to have a dis- 
tinctively weaker interaction with elements in the 'solid'-like phase. Such dif- 
ferentiation is indispensable to the storage of a large amount of energy. Partial 
contact with a heat bath is essential to maintain the differentiation of states 
required to store the absorbed energy. Our coupled pendulum model provides 
a simple example for this behavior. 

The dynamic mechanism for energy storage presented here is simple enough 
to be realized in real protein motors. It is only necessary that the protein 
dynamics exists near a phase transition region with 'solid'-like (vibrational) 
modes and 'gas'-like (rotational) modes and that the folded structure prevents 
some parts of the protein from experiencing the random effect of a heat bath. 
These conditions are expected to be satisfied for many kinds of proteins with a 
sufficiently large size, in addition to the protein motors. It is important to note 
that one can confirm experimentally if the present mechanism of the energy 
storage is valid or not by examining the relationship between the storage time 
and the injected energy, as shown in Fig. ^. 
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Figure 1: Large deviation of the total energy E = K + V at thermal equilib- 
rium, accompanied by the occasional appearance of the rotational mode (see 
the region 50000 <timc< 60000). T = 0.02, AT ^ 10, A^,, = 1 and 7 = lO'^. In 
the upper figure, all time series of the momenta Pi {1 < i < N) are overlaid. 



Figure 2: Process of energy storage after an instantaneous kick at t = to. 
T = 0.02, 7 = 10-2, Nh ^ 1, N ^ 10 and Eo = 0.35. upper: Time series 
of the momenta Pi {1 < i < N). The kicked pendulum A relaxes to thermal 
equilibrium around t^. lower: The time series of the total energy E. The 
dotted line indicates the value at thermal equilibrium. 



Figure 3: Distributions of the lifetime log]^o(^R ~ ^0) for energy storage at 
T = 0.02. A^ = 10 and Eo = 0.5. Each distribution is obtained from 200 
samples. Unshaded distribution: Nh — 1 with 7 = lO^^. Shaded distribution: 
A^,, = 10 with 7 = 10-3. 



Figure 4: Average lifetime {Ir — to) versus the square root of injected energy 
Eo for the system with N — 10. The ensemble average is computed from 100 
samples of different initial conditions, o : Hamiltonian system with 2{K) ~ 
0.02. X : dissipative system with A^^^ = 1, 7 = 10-^ and T = 0.02. A 
: dissipative system with = 1, j = 10~^ and T = 0.02. The data from 
Nh — 9 with the same F value agree with the case N^ — l, within the statistical 
error. 
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N.Nakagawa and K.Kaneko Figure 1 . 




N.Nakagawa and K.Kaneko Figure 2. 




N.Nakagawa and K.Kaneko Figure 3. 
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N.Nakagawa and K.Kaneko Figure 4. 



